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ABSTRACT 

We study theoretical interpretations of the ~ 150-d (superorbital) modulation observed in X- 
ray and radio emission of Cyg X-1 in the framework of models connecting this phenomenon 
to precession. Precession changes the orientation of the emission source (either disc or jet) 
relative to the observer This leads to emission modulation due to an anisotropic emission 
pattern of the source or orientation-dependent amount of absorbing medium along the line of 
sight or both. We consider, in particular, anisotropy patterns of blackbody-type emission, ther- 
mal Comptonization in slab geometry, jet/outflow beaming, and absorption in a coronal-type 
medium above the disc. We then fit these models to the data from the RXTE All Sky Monitor, 
CGRO BATSE, and the Ryle and Green Bank radio telescopes, and find relatively small best- 
fit angles between the precession and orbital planes, 10°-20°. The thermal Comptonization 
model for the X-ray emission explains well the observed decrease of the variability amplitude 
from 1 to 300 keV as a result of a reduced anisotropy of the emission due to multiple scat- 
terings. Our modeling also yield the jet bulk velocity of ~ (0.3-0.5)c, which is in agreement 
with the previous constraint from the lack of an observed counterjet and lack of short-term 
X-ray/radio correlations. 

Key words: accretion, accretion discs - radiation mechanisms: thermal -radio continuum: 
stars - stars: individual: Cyg X-1 - X-rays: binaries - X-rays: stars. 



1 INTRODUCTION 

Periodic variability of emission from Cyg X-1 flux at vari- 
ous frequencies at the period of ~150 d has been reported 
by, e. g., Brocksopp et al. ( 1999a). IPoolev. Fe nder & Brocksopp 
(ll999h.lOzdemir & Demircark (2001). Benlloch et al. t2001..,2004l) . 
iKaritskava et alj i200lh and (Lachowicz et al. ( 200^ hereafter 
L06). This period is much longer than the 5.6-d orbital pe- 
riod terocksopp eTZI IT999h) . and this type of periodicity (or 
quasi-periodicity) in binaries is called superorbital. The gener- 
ally accepted interpretation of the underlying cause of super- 
orbital periodicity in X-ray binaries^ is prece s sion of the ac- 
cretion disc and/or je t (e.g. [X atz 1973, 198^; 'Larwood 1998; 
Wiiers&Pringi^l 19991 ; lOgilvie & Dubus.,2001. ; .Torres et al..,2005. ; 



Caproni et aU |2006|; L06; but with the exc eption of 4U 1820- 



303. e.g.JZdziarsk i. Wen & Gierlinskill2007l) . However, the ques- 
tion arises in which way the precession causes the modulation of 
the observed flux. There appears to be a number of possibilities. 

Considering the X-ray modulation first, the outer edge of the 
optically thick disc may partially cover the X-ray source. This, 



however, would require extreme fine-tuning. Namely, the X-ray 
source has the size ~ lO^iig (where Rg = GM/c?), as indicated 
by the X-ray power spectrum and agreement with theoretical pre- 
diction on the range of radii where most of the accretion power is 
released, while the disc size is generally much larger, up to the or- 
der of the size of the Roche lobe (~ lO'^i^g in Cyg X-1). Another 
possibility is that the outer part of the disc fully obscures the X-ray 
source, but we see the X-rays scattered in a lar ge corona abo ve the 
disc (this appears to take place, e.g., in HerX-l jLe^vl2002[) . This, 
however, would dramatically affect the X-ray power spectrum, re- 
moving oscillations at all frequencies above 1 Hz, which effect is 
clearly not seen, and thus this scenario can be ruled out. The bound- 
free absorption in a spatially extended medium of moderate optical 
depth associated with the outer regions of the disc appears to be 
ruled out as there are a rather weak or no energy dependencies of 
the modulation, see L06 and Section[3]below. (Bound-free absorp- 
tion in the wind from the companion is responsible for th e orbital 
modulation of the X-rays in Cyg X-l. lWen et alj|l999l) . On the 
other hand, a viable scenario is the wind/corona around the outer 
disc being almost fully ionized, with scattering away from the line 
of sight being responsible for the X-ray superorbital modulation. 
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Yet another possibility is that the X-ray emission is in- 
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trinsically aniso t ropic. Such a possibility was considered by 
iBrocksoDP et"ZI ( Il999bh who have relied on the blackbody-type 
anisotropy, where the flux is proportional to the projected area. 
However, there is overwhelming evidence that the dominant ra- 
diative process producing X-rays in the hard state of Cyg X-1 
(and other black-hole binaries) is thermal Comptonization (e.g. 
|Poutanenll998l : IZdziarski & Gierlinskil 20041) . Therefore, it is of in- 
terest to study models of the anisotropy of the Comptonized emis- 
sion to see whether they can reproduce the observed superorbital 
variability. 

We note that the hot Comptonizing plasma most likely also 
forms the base of the jet, which is present in the hard state 
and radiates, at larger distances from the black hole, nonther- 
mal synchrotron radio emissio n correlated with the X-rays (e.g. 
iGallo, Fender & Poolev 1 120031) . However, the radiative process 
giving rise to the X-rays is still thermal Comptonization. Early 
models accounting for the radio-X-ray correlation postulated that 
the observed X-ray emission is due to nonthermal synchrotron 
emission of very energetic power-law electrons with a fine-tuned 
high-energy cutoff. However, that process cannot account, e.g ., 
for the sharpness of the observed cutoffs jZdziarski et al.ll2003h . 
Then, recent X-ray jet models turned to thermal Comptoniza- 
tion to account for the high-energy cutoff. However, the elec- 
tron temperature in that model is very high, kT ^ 3^ MeV 
( iMarkoff. Nowak & Wilmll2005l) . yet still fitted to the cutoff ob- 
served at ~100 keV. Those authors do not explain how it is 
done; if it is due to the 1st order scattering by the thermal elec- 
trons, very strong fine tuning in all hard-states of black-hole bina- 



ries is obviously required. Generally, assuming the jet X-ray ori- 
gin also l eads to a number of other conflicts with the observa- 
tions (e.g.lPoutanen & Zdziarskil2003l : IZdziarski et alj2003ll2004 
iMaccarondboOSl i lYuan etai]|2007h . 



The synchrotron, radio, emission of the jet may be isotropic 
in the comoving frame (in the presence of a tangled magnetic 
field). However, if the jet bulk motion is at least mildly relativis- 
tic, the observed flux will depend on th e jet angle, and will thus 
change with the jet precession. Indeed, [Stirling et alj i200lh find 
that the absence of an observable counter jet requires such a veloc- 
ity for the extended part of the jet. This constraint combined with 
one from the lack of r adio-X-ray correlations on short time scales 
jGleissner etai]|2004h leads to the estimate on the jet velocity of 
~ (0.5-0.7)c. If the jet precesses together with the disc, the ra- 
dio emission will be modulated with the precession period. On the 
other hand, we have no information on the velocity in the core of the 
jet, which can be in principle much lower. In that case, the process 
responsible for the radio superorbital modulation may be free-free 
absorption in the wind from the companion star, which will depend 
on the direction the jet is inclined. This may provide an alternative 
explanation of the superor bital modulation of the radio emission 
dSzostek & Zdziarskill2007l) . Furthermore, both the Doppler beam- 
ing and precession-dependent absorption may take place in Cyg 
X-1. 

Here, we study the precession physical scenarios leading to 
superorbital modulations in a systematic way. In general, any emis- 
sion anisotropic (in the rest frame of the system) with respect to the 
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disc/jet axis will be observed to be modulated when the direction 
of that axis changes. 



2 THE DATA 

We use the X-ray dwell data from the All-Sky Mon- 
itor (ASM) aboard Rossi X-ray Timing Expl orer (RXTE\ 
iBradt. Rothschild & Swan3ll993l ; llxvine et alj|l996l) . We also use 
Compton Gamma Ray Observatory BATSE data in the 20-100 keV 
and 100-300 keV energy ranges, the 15-GHz radio data from the 
Ryle Telescope of the Mullard Radio Astronomy Observatory, and 
the 2.25 and 8.30 GHz data from the Green Bank Interferometer 
(GBI) of the National Radio Astronomy Observatory, Green Bank, 
WV, USA. See L06 for a detailed description of those data sets. 
Compared to the analysis of L06, we also include more recent ASM 
and Ryle data. We study the RXTE/ASM data for MJD 50087- 
53789 (1995 May 1-2006 February 23), the BATSE data for MJD 
48371-51686 (1991 April 25-2000 May 22), the 15 GHzRyledata 
for MJD 50377-53791 (1996 October 10-2006 February 25), and 
the GBI data for MJD 50409-5 1 823 (1996 November 22-2000 Oc- 
tober 6). Hereafter, we refer to the RXTE/ASM channels of 1.5-3, 
3-5 and 5-12 keV as the ASM A, B and C, respectively, and the 
BATSE 20-100 keV and 100-300 keV data as the BATSE A and 
B, respectively. 

Cyg X-I is a highly variable source. Therefore, in order to 
accurately analyze its superorbital variability it is preferable to 
use observations affected in the least way by its aperiodic vari- 
ability. Thus, we study only the hard state, in whi ch Cyg X-I is 
for m ajority of the time. We define it following iZdziarski et al.l 
requiring the average photon spectral index derived from 
the RXTE/ASM fluxes to be < 2.1. In addition, we exclude hard- 
state intervals with high X-ray variability, following the criterion 
used by L06. Namely, we define MJD 50660-50990 as our ref- 
erence interval. We then include only those 30-d intervals of the 
ASM data where < 40 per cent of points exceed by 4(t the av- 
erage flux in the reference interval. This has resulted in removal of 
the following time intervals: MJD < 50350, 50590-50660, 50995- 
51025, 5I400-5I640, 5I840-5I960, 52100-52565, 52770-52880, 
52975-53115, 53174-53554, > 53690. We apply the above selec- 
tion to RXTE/ASM and Ryle data sets, as shown in Fig.[T] 

On the other hand, the GBI data cover a relatively short time 
interval, also shown in Fig. [T] That interval is entirely in the hard 
state. Therefore, we have not applied any additional screening cri- 
teria to those data except the removal of the data marked as bad. 
Then, for the BATSE data, we exclude the periods of the soft state 
of 1994 and 1996, i.e., MJD 49250-49440 and MJD 50230-50307 
(as in L06). 



3 PERIODIC MODULATIONS IN THE DATA 

The currently most comprehensive analysis of the orbital and super- 
orbital modulations of Cyg X- 1 has been done by L06. In particular, 
they found the superorbital period of ~150 d consistent with con- 
stant in all the available data since ~1976. Furthermore, they found 
the phase of the superorbital also compatible with constant (see fig. 

7 and table 3 in L06). 

Here, we use the Lomb-Scargle method jLomblll976l : IScarglel 




Figure 2. The Lomb-Scargle periodograms for (a) the ASM C and (b) 
the Ryle data (before prewhitening; corresponding periodograms after 
prewhitening are shown by L06). 



both data sets. In addition, we see strong peaks corresponding to 
periods longer than 150 d for the Ryle data. Their origin remains 
unclear though they are likely to be artefacts, see L06 for discus- 
sion. 

We first consider the orbital modulation (caused by absorption 
in the wind. IWen et al.lll999l : iBrocksopp. Fender & Poolevll2002l: 
Szostek & Zdziarski '20071; see the latter paper for corrections to 
jrocksopp et al. 2002). For it, we use the sp ectroscopic ephemeris 
terocksopp et al.lll999al : [LaSaIa et al.lll998 



min[MJD] = 50234.79 + 5.599829S, 



(1) 



where E is an integer. At those times, the companion star is in 
front of the X-ray source. We then fit the light curves folded over 
the orbital period and averaged within each of 20 phase bins. We 
use then logarithm, G = In F, of the photon or energy fluxes or 
count rates, F, for fitting with the sum of harmonics, 



Gn.od(<^) = (G) - ^ Gk cos[27rfc( 



(2) 



1 1982 ) to characterize the periodicities in the data analyzed by us. 
The periodograms for the ASM C and Ryle data are shown on Fig. 
|2] We clearly see the peaks corresponding to 5.6 d and ~150 d for 



where (p is the 0-1 phase, (G) is the (fitted) average value of the 
logarithm of the flux, and Gk and <l)k are the amplitude and the off- 
set phase, respectively, for the fc-th harmonic. We adopt the conven- 
tion that Gk > and —0.5 < ktfjk < 0.5. The obtained parameters 
for Ai' = 3 for the RXTE/ASM and Ryle data and iV = 1 for the 
GBI data (see L06) are given in Table [T] We do not consider here 
the orbital modulation for the BATSE data since it is very weak, 
< 3 per cent (L06). 

We note that assigning errors to the folded and averaged fluxes 
is not a unique procedure. Here, we first divide the light curve 
into time bins of the 1/20 of the period (orbital or superorbital) 



© 2007 RAS, MNRAS OQO.fTlflOl 



4 A. Ibragimov, A. A. Zdziarski and J. Poutanen 

Table 1. The coefficients for the orbital modulation fitted by equation (2)- The units of exp(G} are s^^ for the RXTEI ASM data, and Jy for the Ryle and GBI 
data. 



Detector 


(G) 






G2 


02 




</.3 


RXTEI ASM A 


1.829±0.005 


0.113±0.008 


0.007±0.011 


0.042±0.008 


-0.02±0.01 


0.017±0.008 


-0.01±0.02 


RXTEI ASM B 


1.825±0.004 


0.066±0.006 


0.005±0.013 


0.022±0.006 


0.03±0.02 


0.008±0.006 


-0.00±0.04 


RXTEI ASM C 


2.188±0.003 


0.035±0.005 


0.005±0.022 


0.011±0.005 


-0.04±0.04 


0.007±0.005 


-0.02±0.04 


Ryle 


^.429±0.007 


0.167±0.011 


0.14±0.01 


0.035±0.010 


0.16±0.02 


0.017±0.010 


0.15±0.03 


GBI (8.30 GHz) 


^.122±0.011 


0.064±0.015 


0.17±0.04 










GBI (2.25 GHz) 


^.187±0.010 


0.024±0.014 


0.4 ±0.1 











Table 2. The periods and amplitudes of the superorbital modulation in the analyzed data. The amplitudes have been calculated using the average superorbital 
period of L06, 151.43 d. The errors are Icr. The second set of the A values gives the modulation amplitudes obtained by L06. 



Detector 


Energy 


P[d] 






A 




A (L06) 


RXTEI ASM A 


1.5-3 keV 


153.5 


± 


1.5 


0.120 


± 


0.005 


0.133 


± 


0.005 


RXTEIASM B 


3-5 keV 


153.5 


± 


1.5 


0.105 


± 


0.003 


0.114 


± 


0.003 


RXTEIASM C 


5-12 keV 


153.5 


± 


1.3 


0.103 


± 


0.003 


0.114 


± 


0.003 


BATSE A 


20-100 keV 


151.0 


± 


1.1 


0.083 


± 


0.001 


0.081 


± 


0.001 


BATSE B 


100-300 keV 


152.0 


± 


3.0 


0.075 


± 


0.002 


0.070 


± 


0.001 


Ryle 


15 GHz 


152.7 


± 


1.4 


0.140 


± 


0.005 


0.105 


± 


0.002 


GBI 


8.30 GHz 


148.8 


± 


3.1 


0.116 


± 


0.004 


0.123 


± 


0.003 


GBI 


2.25 GHz 


150.7 


± 


3.5 


0.122 


± 


0.004 


0.107 


± 


0.003 



and obtain the average of the flux, Fij, where i is the number 
of the phase bin and j is the number of the time bin contribut- 
ing to the i-th phase bin, weighted by the inverse squares of their 
measurement error. In this way, we avoid any contribution to our 
folded/averaged light curves from the source variability on time 
scales shorter than that corresponding to the length of our cho- 
sen phase bin. The resulting measurement uncertainties on Fij are 
rather small, and we neglect them. Then, we consider the aver- 
age and standard deviation of all time bins contributing to a given 
phase bin. Here, we want to take into account the actual aperiodic 
variability of the source, not just the measurement errors. There- 
fore, we calculate the unweighted average within the i-th phase bin. 



Fi = (^X]j=i Pijj /-^i' where Ni is the number of time bins con- 
tributing to the i-th phase bin, and the corresponding rms standard 
deviation, a„ given by af = [Ef^i(^'« " F,jf]/[N,{N, - 1)]. 
Then, we take Gi = in Fi and fit with equation 

Subsequently, we prewhiten the logarithmic light curves by 
subtracting from them the orbital variability fitted above. 



G'(f) = G(t)-Gn.odW + (G). 



(3) 



We now again apply the Lomb-Scargle method to the prewhitened 
data. (We do not prewhiten the BATSE data, as their orbital mod- 
ulation is very small, see above.) The obtained main (~150-d) su- 
perorbital periods are presented in Table |2] The given la errors 
have been estimated us i ng the method described in section 4.4 of 
ISchwarzenberg-Czemvl jl99lh . The superorbital period of Cyg X- 
1 was determined by L06 using the weighted average, 151.43 d, of 
periods determined for a large number of data sets. For our data 
alone we find 152.4 ± 0.6 d, and including the Ginga and Ariel 5 
data from L06, we obtain a similar value, P — 152.0 ± 0.4 d (la 
error), consistent with that of L06. Thus, for consistency with L06, 
we hereafter use their ephemeris, 

min[MJD] = 50514.59 + 151. 43S. (4) 
We then fold and average within each of 20 superorbital phase 
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Figure 3. The superorbital phase diagrams for two of our data sets, showing 
the light curves folded over the superorbital ephemeris, equation @> ^nd 
averaged within each phase bin. (a) The ASM A X-ray data modelled by 
thermal Comptonization in the slab geometry. The best fit coiTesponds to 
i = 32° and 5 = 18° (see Section|42] Fig.|9). The Ryle radio data fitted 
with the jet model (see Section l4~n and Fig.|6), with the assumed values of 
i = 40° and /3 = 0.4, which yield the best-fit of 5 = 11° . 
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observer 




Figure 4. The assumed geometiy of Cyg X-1 at the zero orbital phase (i.e., 
at the superiour conjunction of the X-ray source). 

bins the prewhitened data, and tit them using equation |2) with A'^ — 
1. The errors used in the fit are calculated in the way described 
above. The fits give us the amphtude of the variability, which we 
define as, 

^ Fmax ^min R 1 

Jmax + i*min i? + 1 ' 

where i^max and Fmin are, respectively, the maximum and min- 
imum fitted fluxes, which ratio is 7? = Fmax/^min, and R = 
(1 + A)/{1 — A). For our harmonic fit with A'^ = 1, Fmax.min ~ 
exp((G) ± Gi). The obtained values of A are hsted in Table|2] 
We see that the amplitudes of the superorbital modulation in both 
the ASM and radio data are almost independent of the energy, con- 
firming the previous result of L06. On the other hand, we confirm 
the result of L06 that the modulation amplitudes in the 20-300 keV 
energy range are significantly smaller than those in the 1.5-12 keV 
range. Two of the resulting superorbital phase diagrams are shown 
in Fig. [3] 

We point out that the obtained amplitudes are somewhat af- 
fected by systematic uncertainties, in particular those related to the 
selection criteria. We have tested that indeed changing the criteria 
used to choose the fitted parts of the light curves results in some 
changes of the values of A. To illustrate this effect. Table [2] also 
gives the values of A of L06, who used both somewhat differ- 
ent selection criteria and different prewhitening method (namely, 
prewhitening the data with all other significant modulations, not 
only the orbital one). We can see that their values differ from ours 
by typically 10 per cent, which can be taken as the systematic un- 
certainty of the values of A. The largest difference is for the Ryle 
data, which appears to result from not applying the selection crite- 
ria obtained from the ASM data to the Ryle data in L06. Also, we 
note that L06 treated each modulation as due to absorption, which 
lead to increases of the average fluxes after each prewhitening. In 
our method, we preserve the average fluxes, see equation ([3}. 



4 FLUX VARIATION FROM PRECESSION 

Let us consider a precessing slab model with a constant tilt angle, 
S, see Fig.|4l and a constant angular velocity lo = where 
P is the superorbital period. Then, the superorbital phase is propor- 
tional to time within each cycle, e.g., (pso ~ t/P for t < P. Note 



that the precession is independent of the orbital motion of the sys- 
tem, i.e., the disc inclination in an inertial frame does not directly 
depend on the orbital phase. The precession may also include any 
structure attached to the disc, e.g., a corona and/or a jet. The an- 
gle, tf), between the normal to the disc and the direction towards the 
observer is given by 

cos tp — I cos S cos i — sin 5 sin i cos(27r(?!>ao)|, (6) 

where i is the orbital inclination to the line of sight. Here, (j>so = 

corresponds to the disc being most inclined with respect to the 
observer. The emissivity of the disc (and any additional structure, 
e.g., a jet) integrated over its surface and over the orbital period can 
be only a function of the direction with respect to the normal, i.e., 
tp. The value of varies between i + S at ipso ~ and |i — (S| at 
(pso = 0.5. We note that since no secondary maxima are observed 
in the X-ray folded superorbital light curves (Fig.jSj, we see only 
one side of the disc, i.e., always tp < 90°, implying i + S < 90°. 
Otherwise, we would see at some moment only the disc edge and 
then its opposite side. Note also that equation (|6} is symmetric with 
respect to an exchange of i and S. Thus, for any configuration with 

1 > 5 there is an equivalent one with 5 > i but with the two values 
exchanged. Thus, we will consider only the case of i > 5, keeping 
in mind the existence of the symmetry. 

4.1 Analytical models 

If the emission from the precessing structure is anisotropic, flux 
variations over the precession period take place. We consider first a 
few simple models. Here, 0° < ip < 90° . 

(a) F{tp) = Gcos^. This case corresponds to a slab/disc 
emitting with a constant specific intensity, i.e., blackbody-like. 

(b) F{'tp) — Ccosip{l + acostp). A more general case, 
where a is an anisotropy parameter parameterizing departures from 
the blackbody. Here, a = corresponds to the constant specific 
intensity (blackbody) case, a = 2 is for Chandra sekhar-Sobolev 
optically-thick elect ron scattering atmosphere iChandrasekhail 
ll960l : ISobolevlll963l) , and some values in the — 1 < a < range 
correspond to Comptonization in a slab with a moderate optical 
depth (Sunvaev & Titarchuk 1985; Viironen & Poutanen 2004'; see 
also Section l4~2l below). Values of a < —0.5 correspond to the 
Thomson optical depth across the slab t < 1. 

(c) Fi'ip) = C[7(l - /3cosi/))]-(^+^\ where F is the en- 
ergy flux. This law corresponds to the emissivity pattern of a steady 
Jet or disc outflow when the retardation effe ct is absent (see e.g. 
iRvbicki & Lightmaril979l ; ISikora et aljl997l) . Here P = v/c is the 
bulk velocity in units ofc,7 = (1 — /9^)~^''^, and F is the observed 
photon spectral index. Note that if this model is applied to the jet 
case, the counter jet may be visible as well (unless obscured by 
the di sc and/or free-free absorbed in the wind. lSzostek & Zdziarskil 
l2007h . In the counter jet is visible, a second term with /3 replaced 
by — /3 should be added to the above expression for F{xp). 

(d) F(V') = Ce~^/°°"''', where r is the optical depth. This 
law corresponds to anisotropy caused by an absorber with the same 
symmetry as the disc. 

Using these models, we relate either the amplitude. A, or the 
flux ratio, R, to the model parameters. We assume 5 < 90° — z, in 
which case the global flux minimum corresponds toip = i + S. For 
the case (a), we have, 

A = tan S tan i. (7) 
In the case (b). 
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Figure 5. The relations between amplitude A, the precession angle, 5, and other parameters of our analytical models for i = 40° . The panels (a-d) correspond 
to the models (a-d) desciibed in Section l4~I1 (a) The dependence (TJ between A and <5. (b-d) The dependencies between model parameters and the precession 
angle <5. The assumed values of the variability amplitude of A = 0.10 and 0.14 are shown by the soUd and dashed curves, respectively. In the case (c), F = 1 
was used. We also show the results at F = 1.7 for the above two values of A by the dotted and dash-dotted curves, respectively. 



A = 



(1 + 2a cos 5 cos i) sin 5 sin i 
cos 5 cos i (1 + a cos 5 cos i) + a sin^ 5 siv? i ' 



and 



sin S sin i — A cos S cos i 
Acos'^ 5 cos^ i — 2 cos S cos i sin 5 sin i + Asin'^ 5 sin^ 



(8) 



.(9) 



Note that in equations we have assumed that the flux 

monotonously decreases with increasing disc viewing angle, 
This is generally the case for a > —0.5. If a < —0.5, the above 
formulae apply only in the range of > arccos( — l/2a). Other- 
wise, the maximum flux is achieved either altp — arccos(— l/2a) 
or at the maximum value of -i/; of the precession (i.e., i + 5). In the 
former case, there will be two local minima of the flux, and in the 
latter, the minimum flux would be achieved at the minimum value 
of ^. Both those cases can be ruled out for Cyg X-1. First, only one 
minimum per superorbital period is seen, and second, the X-ray su- 
perorbital modulation is in phase with the radio one (L06), and it is 
highly unlikely that our case (b) with a < —0.5 would apply to the 
jet radio emission. Therefore, we hereafter assume that that the flux 
monotonously decreases with increasing ^ and equations ([8}(9]( are 
applicable. If anisotropy is interpreted as resulting from Compton 
scattering is a slab, the above constraints on a also mean that the 
slab Thomson optical depth is t ^ 1, which is consistent with the 
estimates of r from the X-ray spectra (see Section l4~2l ). 

In the case (c), we can obtain formulae relating the flux ratio, 
R, to /3 (neglecting hereafter the counter jet emission), 

'1 -/3cos(i + ^)"i^+^ 



R = 



- /3cos(i — S) 
j^i/(i+r) 



i?i/(i+r) cos(j - 5) - cos(i + 5) 



(10) 



(11) 



(see also lGregorv et aLlligSg *). Note that for given R and i, there is 
a low limit on 5 (as /3 < 1). The minimum /3 required to explain 
an observed flux ratio at a given inclination is achieved when the 
precession angle is given by 



i?i/(i+r) + l 
t^"^^ j^i/(i+r) „i ^'^i'- 



(12) 



Note that this value of S is generally > i, i.e., the interchange sym- 
metry between i and 5 is no longer present here. 
In the case (d), we find, 



R = exp {r[sec(i + 5) 
InR 



sec(i — S)]} , 



sec(i + 5) — sec(i — S) 



(13) 
(14) 



Various authors have obt ained different constr ai nts on the 
inclination of Cyg X-1, e.g.. lOies & BoltonI jl986h, IWen et all 
(1999), Abub ekerov. Antokhina & Cherepashchuld j2004l) and 
IZi61kowskil ( l200^ goti = 28°-39M0°-40°, 31°-44°, and 23°- 
38°, respectively. Hereafter in this Section, we assume i = 40°. 
Relationships between S and other model parameters are shown in 
Fig.H] 

We then fit the anisotropy laws (a-d) to the observed profiles 
of the superorbital variability (Section |3]l. The (blackbody-like) 
model (a) yields 5 = 8?2 ± 0?5, 7?2 ± 0?5, 7?0 ± 0?5, 9?5 ± 1?5, 
8? 4 ± 1? 7 and 5 = 8?0 ± 1?6 for the ASM A, B and C bands, 
Ryle and GBI 2.25 and 8.30 GHz data, respectively. For the mod- 
els (b-d), the correlations between the parameters (see Fig. [5) lead 
to elongated error contours, which we show in Fig. |6] The values 
of the statistic xt = X^/dof for the best fits are as follows for the 
ASM A, B and C bands, the 15 GHz Ryle data and the 2.25 and 
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Figure 6. The 90% confidence contours for the models (b-d) (marked on the panels) fitted to the observed superorbital variability profiles (Section[3) at the 
assumed i = 40°. In the case (c), we use F = 1.7 and 1 for the fits to the X-ray and radio data, respectively. The soUd and dotted contours are for the ASM 
A data and the Ryle data, respectively. The ASM B, C and GBI data yield similar respective contours. 



8.30 GHz GBI data, respectively. For the model (a): xt = (8/17, 
11/17, 14/17, 10/17, 11/17, 21/17), (b): (8/16, 11/16, 14/16, 10/16, 
11/16, 17/16), (c): (9/16, 11/16, 13/16, 11/16, 12/16, 22/16), (d): 
(8/16, 11/16, 14/16, 10/16, 11/16, 17/16). Thus, we see that we 
cannot distinguish between the models based on the fit quality. 

For the model (c), we use F = 1.7 (the average X-ray power- 
law index of the hard state of Cyg X-1) for the ASM data, and 
r = 1 (corresponding to the obs erved 2.2-22 GHz radio emis- 
sion of Cyg X-1 in the hard state. iFender et al. 2000) for the Ryle 
data. In the latter case, we do not include emission of the counter 
jet. The results for those two cases are consistent with a coro- 
nal outfl ow at /3 c:^ 0.3 explaining the X-ray Co mpton reflection 
strength (Malz ac. Beloborodov & Poutanenl200 l'). and the inferred 
radio jet velo city of /3 ~ 0.5-0.7 in Cyg X-1 (istirling et al. 2001 ; 
ICjleissner et alj2004l) . respectively. The model (c) fitted to the Ryle 
data superorbital phase diagram is shown in Fig.[3tb). 

We point out that the jet in Cyg X-1 is embedded in the stel- 
lar wind from the companion, which causes the orbital modulation 
of the radio emission via free-free absorption ( Szostek & Zdziarskil 
l2007h . Thus, a jet precession, which changes the path of the jet dur- 
ing the course of the orbital motion, will also change the average 
optical depth in the wind along the line of sight. This will give rise 
to a superorbital modulation in addition to that caused by the jet 
beaming. 



4.2 Anisotropy of thermal Comptonization 

The dominant radiative process giving rise to X-rays in the 
hard state of Cyg X-1 is, most likely, Comptonization of some 
soft seed photons by predomi nantly thermal electrons at the 
temperature of ~50-100 keV jGierliriski et alj 1 19971 : IPoutanenI 



199S 


; Di Salvo etal 2001: 


Frontera et al 2001 ; Zdziarski et al.l 


2002 


; iMcConnell et al.ll2002 


; llbrasimov et al.ll2005l). This is evi- 



denced mostly by the characteristic form of its high-energy cut- 
off, present also in oth e r black-hole binaries in the hard state 
(e.g. iGrove et al.' '1998"; 'Zdziarski et alj 1 19981 : IWardziiiski etall 
I2OO2I : .Zdziarsk i & Gierlinski 2004) . The location of the ther- 
mall y Comptonizing plasma is eithe r a hot inner accretion flow 
(e.g. 'Poutanen, Krolik & Rvde' '1997"; 'Esi n et alj 1 19981 : iPoutaned 
1998; Zdziarski et al. 2002; Yuan et al. 20o3) or a coronal outflow 
( lBeloborodovlll999l ; lMalzac et al.ll200ll) . 

Thus, we fit such a model to a representative hard-state 
spectrum of Cyg X-1. We choose the RXTE observation 10238- 
01-03-00 toget her with the CGRO/O SSE observation VP 612.5 
(spectrum 6 of llbragimov et all |2005| ). which gives us a broad- 
band energy coverage of ~3-1000 keV. Our detailed model 
consists of the main component due to thermal Comptoniza- 
tion in a sla b geometry calculated using the iterative scatter- 
ing method of IPoutanen & SvenssonI ( 119960 (COMPPS in XSPEC, 
Arnaud '199^. In addition, we include Compton reflection 
( Magdziarz & Zdziarski 1995), a fluorescent Fe Ka line, and a 
soft excess. We attribute the last component to thermal Comp- 
tonization as well, but in a plasma with the Compton y parameter 
much smaller than that of the main C omptonizing plasma, and ne- 
glecting its Compton reflection (as in lFrontera et al.|[200ll) . Here, 
y = ArkTc/nicf? , where r is the optical depth of the plasma, and 
Te is its electron temperature. The seed photons for Comptoniza- 
tion are from a blackbody spectral component at the temperature of 
Tbb, which we keep equal to 0.2 keV in the fit. Also, we assume the 
inclination of = 40° . The spectrum with the resulting fit compo- 
nents is shown on Fig.|7] The main fit parameters are fcTc = 88 keV 
and r = 1.17 for the main Comptonization component, the solid 
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Figure 7. An example of the broad-band X-ray spectrum of Cyg X-1 in 
the hard state (data points), see Section 14.21 The data are fitted by a ther- 
mal Comptonization and reflection model. The solid, dashed, dot-dashed 
and dotted curves represent the total model, the main Comptonization com- 
ponent, the soft excess (also fitted by Comptonization), and the reflection 
component together with an Fe Kq fluorescence line, respectively. At low 
energies, the spectrum is absorbed by interstellar and circumstellar media. 
The bottom panel shows the residuals to the fit. 



gure 9. 1 he 90 per cent confidence contours tor the inclination and 
precession angle. The solid and dot-dashed curves correspond to the joint fit 
of all thi-ee ASM data and the BATSE A data, respectively. The dotted and 
dashed curves show the confidence contours of the fit to the Ryle data for the 
jet model with F = 1 for /3 = 0.3 and 0.5, respectively. The two horizontal 
lines mark i = 30° and 45°, approximately corresponding to the range of 
the inclination of Cyg X-1 allowed by other observational constraints. The 
line in the bottom right corner corresponds to i = 5, below which there is 
a symmetric, i < 5, solution (ruled out for Cyg X-1 by the 30 < i < 45° 
constraint). 



2.0 I I I I I I I I I I I I I I I I I I T 




cos 'ijj 

Figure 8. Anisotropy patterns of thermal Comptonization in a slab, obtained 
for our representative X-ray spectrum of Cyg X-1 in the hard state, Fig.|2] 
The sohd, dotted, dashed, dot-dashed, and three-dot-dashed curves give the 
photon flux integrated over the 1.5-3 keV, 3-5 keV, 5-12 keV, 20-100 keV, 
and 100-300 keV energy bands, respectively. The curves are normalized 
to (F) = F(cosV)d cos il>. The blackbody emission pattern would be 
represented by a diagonal straight line, whereas the anisotropy law of the 
case (b) of Section l4~l1 with a = —0.4 lies between the shown dependencies 
for 5-12 keV and 20-100 keV. 



angle of the Compton reflector of n/27r = 0.17, and r = 0.87 
of the soft-excess Comptonizing plasma, for which fcTc = 20 keV 
was assumed. The absorber column density is A'^h ~ 2.3 x lO'^'^ 
cm"^. The fit statistic is good, xl = 371/401. 

We then vary the inclination of this model, and for each ^jj 
we integrate the photon flux in five energy bands, 1.5-3 keV, 3-5 
keV, 5-12 keV, 20-100 keV, 100-300 keV. The resulting depen- 
dencies of F{'4>) are shown in Fig. [8] We also compare the re- 
sulting functions to = cos7/;(l + a cos ?/))/(l/2 + a/3), 
i.e., the normalized case (b) of Section [43] We find a ~ —0.4 
as the closest overall approximation for the shown angular depen- 
dencies. On the other hand, the normalized blackbody dependence 
equals F^t})) — 2 cos tp, which yields less (more) emission than the 
normalized thermal Comptonization dependencies at cos < 0.65 
(>0.65). 

We now use the dependence, Fixjj), in an energy band together 
with 5, (f>Bo) of equation ^ to fit the F{(f>so) observed in that 
energy band, with i and 5 as the fitted parameters. We use here the 
ASM and BATSE data. (Since the BATSE data are given as energy 
fluxes, we use the corresponding integrated energy fluxes for them.) 
This yields very good fits to all the data, e.g., at the best fits to 
the ASM A, B, C data, and the joint fit for all three detectors equal 
9/17, 11/17, 13/17 and 34/57, respectively. We have calculated the 
90 per cent confidence contours for all the data, and found that 
they remain in excellent mutual agreement. We show the contours 
for the ASM joint fit and for the BATSE A data in Fig. [9] by the 
solid and dot-dashed curves, respectively (the BATSE B contours 
are very similar to those of the BATSE A, except for extending to a 
somewhat higher value of 5 at a given i). 

This agreement provides a strong argument for the correctness 
of our model of anisotropy as due to thermal Comptonization in 
a geometry similar to that of a slab. Namely, we see in Table |2] 
that the fractional X-ray modulation decreases monotonically with 
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the increasing photon energy band, from 1.5 keV to 300 keV. In 
particular, the ASM A channel shows the amplitude higher than that 
of the channels B and C, and the BATSE data show the amplitudes 
lower than those of the ASM. This is in excellent agreement with 
the theoretical anisotropy patterns shown in Fig. [8] Namely, we see 
there that the degree of anisotropy of thermal-Compton emission, 
as measured by the slope of the curves, decreases monotonically 
with the increasing energy (for cos i/i > 0.3, which is satisfied in all 
of our fits). This effect is due to the average number of scattering 
in the plasma increasing (thus leading to less anisotropic emission) 
with the photon energy. This then leads to a corresponding decrease 
of the relative amplitude of the precessional modulation, which is 
both present in the X-ray data sets (Table|2}, as well it is very well 
fitted by our model. 

On the other hand, we see the contours are very elongated, 
which results from a strong correlation between i and S in the 
model. The la contours are only slightly smaller than the shown 
90 per cent ones, indicating that the minima of are very shal- 
low. Therefore, in order to constrain 5, we use the constraint of 
30° < i < 45° (approximately corresponding to the constraints on 
i by various authors listed in Section |4Tt . with the boundaries of 
this region shown by the horizontal lines in Fig. [9l Then, we find 
8° <(5< 11° at i = 45°, and 17° <(5<22° at i = 30°. We then 
compare these results with the 90 per cent confidence contours for 
the radio emission as fitted with the jet model. We find that model 
is compatible with the Comptonization model for the X-ray data if 
/3 = 0.3-0.5 (shown by the dotted and dashed contours for /3 = 0.3 
and 0.5, respectively), but not if 13 = 0.7. For (3 = 0.3-0.5 and 
30° < i < 45°, the radio data provide constraints on 5 less restric- 
tive than the X-ray data. We show the best-fit ASM A model and 
a model for the Ryle data on the superorbital phase diagrams, see 
Fig. [3] 

The main caveat to our results here is related to using the 
slab model for Comptonization. It is unlikely that the actual X-ray 
source in Cyg X- 1 has exactly this geometry, and it is probably less 
anisotropic. This would yield somewhat higher precession angles 
than those obtained by us, but it is unlikely to affect our qualita- 
tive conclusions. Also, we have used a specific (but rather typical) 
single spectrum to obtain F{ip). Howe ver, though there is so me 
spectral variability in the hard state (e.g. Ilbragimov et al.|[2005l) , it 
is aperiodic and its effect will be averaged out over a number of su- 
perorbital cycles. Thus, using this assumption is unlikely to affect 
significantly our results. 



5 DISCUSSION AND CONCLUSIONS 

We have studied a general problem of how precession can affect 
the observed fluxes and spectra. Apart from an obvious case of a 
full or partial covering of the central source by an accretion disc 
(e.g., in Her X-1, Leahy 2002), the observed flux can be modulated 
if the emission pattern of the precessing source is anisotropic. We 
have studied a number of likely anisotropy patterns, due to either 
intrinsic anisotropy of the emission pattern in the source rest frame, 
or Doppler boosting of emission isotropic in the rest frame, or ab- 
sorption/scattering by a corona covering the disc. We have derived 
formulae for the flux as a function of the precession phase for those 
processes. In addition, we have numerically calculated the X-ray 
anisotropy pattern due to thermal Comptonization and Compton 
reflection in a plasma with parameters characteristic to black-hole 
binaries in the hard state and assuming a slab geometry. 

We have then applied these results to the superorbital modu- 



lation of the 2.25-15 GHz radio and 1.5-300 keV X-ray emission 
of Cyg X-1 in the hard state. We find that most models make sim- 
ilar predictions for the precession angle, i.e., the tilt between the 
disc plane and the orbital plane, of the order of 5 ~ 10°-20°. The 
statistical quality of all the models is similar, thus, they cannot be 
distinguished on the basis of their x^- The model with absorption 
by a disc corona yields strongly anticorrelated the optical depth and 
the precession angle. 

Our most physical model appears to be thermal Comptoniza- 
tion for the X-rays and jet Doppler boosting for the radio emis- 
sion. The joint constraint to the X-ray and radio data yields for 
/3 = 0.3 the values of 5 > 10°, «<47°; and for /3 = 0.5, the 
values of 5° < (5 < 10°, 40° < i < 60° (see Fig. ^. If we assume 
30° < i < 45° (from other observational constraints on the in- 
clination, see Section l4Tt . we obtain 8° < 5 < 22° from the X-ray 
data, which then implies (3 ~ 0.3-0.5 using the Ryle data. Thus, 
both the disc/jet precession angle and the jet velocity in Cyg X-1 
appear relatively small. We also point out that our Comptonization 
model also explains and well fits the decrease of the fractional su- 
perorbital modulation amplitude with the increasing photon energy 
observed in the 1.5-300 keV X-ray emission of Cyg X-1. The de- 
crease is explained by the Comptonization anisotropy decreasing 
with the photon energy due to the correspondingly increasing aver- 
age number of scatterings suffered by a photon before leaving the 
source. 

Finally, it is interesting to compare our results for Cyg X-1 
with those for the well-known precessing jet source SS 433. In that 
case, the jet velocity and the precession angle are very precisely de- 
termined based on the Doppler shift of spectral li nes (formed in the 
jets). Thei r values are P ~ 0.265 and S ~ 21° jEikenberrv et al.l 
I2OOII : see lFabrikall2004l for a review). They are rather similar to 
those determined by us for Cyg X-1. However, the disc emission 
is not seen at all due to heavy obscuration, and the radio emission 
shows no periodic variability on either orbital or precesssional pe- 
riod. This appears to be related to the rather large inclination of SS 
433 of j ~ 78°, at which the disc/jet system goes through the edge- 
on orientation at some point of each precessional cycle. Then, the 
counter jet is seen at a lower inclination than the primary jet dur- 
ing a part of the cycle, and both of them contribute comparably 
to the observed flux (see Section |4Tt . Also, free-free ab sorption 
of the radio emission is of importance in that system (e.g. lFabrikal 
'2004*), which introduces a major complication to predictions of the 
precessional variability given that details of the wind geometry are 
unknown. On the other hand, the precess ional modulation is very 
significant in the X-rays jWen et al.ll2006l) . emitted by the jets. 

The disc/jet precession in SS 433 is probably caused by the 
precession of the companion star, with the rotation axis inclined 
with respect to the orbital plane (see Fabiika 2004 and references 
therein). On the other hand, this precession mech anism is unlikely 
in Cyg X-1, where the orbit is strictly circular jBrocksopp et al.l 
Il999b [). indicating that tidal forces have also brought the stellar ro- 
tation to alignment. Then, the precession may be caused, e.g., by 
the gravitation al force of the star acting on the inclined disc (e.g., 
iLarwoodl 19981 : L06). 
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